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irrRC DUCT ION 


The Nttwton-Kiintorovich technique ia a rijjorouJ nuthenut leal concept 

baseil on functional nnalywis that transfornu u non-linear partial <llf- 

fercnLlal equation into a sequence nf linear partial differential 

equations whose solutions converge to the solution of the original 

non-linear problem providing appropriate conditions can be satisfied. 

It has been employed operationally by Bellman^ to solve a variety of 

non-linear problems with considerable success, and vms more recently 

2 

used by Morihara and Cheng to produce a numerical solution of the 

3 4 

steady-state Navler-Stokes equations in two-dimensions. Gabrielsen ' 
has established a mathenatical justification for the use of the tech- 
nique in solving the non-steady Navier-Stokeu equations in two dimensions. 

The Newton-Rantorovich technique has major advantages erver other 
techniques used to solve the non-steady Navier-Stukes equations. The 
main computational advantage is t>iat a sequence of linear partial 
differential equations need only be solved rather than a single non- 
linear one. A linear partial differential equation can be solved 
numerically by simply solving a system of linear algebraic equations 
directly, whereas a non-linear equation requires an elaborate iteration 
procedure that may require more compute time and wiiosc convergence is 
often difficult to achieve especially at high Reynolds numbers. A 
theoretical advantage of the technique is that explicit error estimates 
are attainable (that is, the difference between each linear solution 
and the solution to th' original non-linear problem). 
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The obvious disadvantage of the technique is tliat more than one 

linear partial differential equation must he solved li order Cr achieve 

a solution. In fact, based on the results of previous Investigations 

of the technique, the number of iterations of the procedure required 

for convergence of the technqiue would appear to be verv problem depen- 

2 

dent. Moriharn and Cheng applied the technique to a steady flow in 
the entrance region of a straight channel for various Reynolds numbers 
and found that *'.e nuniber of iterationr required for convergence increased 

3 

rapidly with Reynolds number. In addition, Gabrielsen established an 
upper bound for the error associated with the technique as applied to 
the non-steady Navler-Stokes equations as a function of the number of 
successive iterations and the Reynolds number. His analysis indicated 
that the number of iterations required for convergence depended on the 
accuracy of the initial guess at the solution, and that an Increasingly 
more accurate initial guess at the solution is required with Increasing 
Reynolds number in order to achieve convergence of the procedure in the 
same number of Iteraticns as required at lover Reynolds numbers. This 
Reynolds number dependence is due to the fact that the solution to the 
Navler-Stokes equations becomes progressively more non-uniform and 
more concentrated in a smaller region of the flow field with increasing 
Reynolds number. The convergence of the procedure must also depend on 
the shape of the geometry under consideration because the severity and 
location of these non-unifoinnl ties and concentrations will vary with 


geometry. 


For expediency, the actual Impleitentatlon of the procedure tms 
acconplishcd by modifying an exiatlng computer program vrrltten by 
Hehta^. Mehta *h code 1 h considered one of the cost advanced computer 
programs in existence for numerically eulving Che laminar, unsteady. 
Incompressible davier-Stok.es equations in evo dimensions, in that 
the program is capable of determining the flow field about a flat 
plate, a circular or elliptical cylinder, and synnetric or cambered 
airfoils at arbitrary angles of attack including stall. The code is 
also capable of performing solutions at ary Reynolds nmber for which 
the laminar flow assumption is reasonable. Another feature of the 
prograci that was particularly ureful for Cite current study was that 
the equations of motion were transformed into the Interior cf the unit 
circle and the coordinate perpendicular to the surface of the geometric 
shape was stretched depending on the Reynolds number so chat the solution 
produced by the transformed equation was relatively similar for all 
geometries and Reynolds numbers. This feature also produces a trans- 
formed Navier-Stokes equation whose solution is relatively uniform 
througitout the region of calculation. 

For the purposes of our investigation, Mehta's code was modified 
to solve each linear partial differential equation produced by the 
Newton-Kantorcvich process. In ad<Iitlon, due to computer storage 
limitations, the program had to be simplified to perform computations 
for only symmetric geometries and flow fields. Test cases were performed 
for a circular cylinder at a Reynolds number of 15 and *or a syninetrlcal 
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122 thick airfoil at zero angle of attack at Keynnjcla nutobera of 10^, 

4 5 

10 , and 10 . The rcaulta of the calrulatlona ualni; the Newton-Kantoiovlch 
procedure were In all canes coopared with the results obtained foi: the 
same problem by the original unmodified version of the program. 
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MATHEMATI'JAL FORMULATION 


Th« flow field exterior to the deoa^ttv under coneideretion is 
■sapped into the unit circle. As a result of this transformation, 
the equation governing the unsteady, incompressible flow of a Newtonian 

fluid may be expressed in terms of the vorticity w and the stream function 

u. 7 
Y as 

- i ^ I'r' ^ t ^ r‘) - * ^ - r if ^ ) a) 

"'Li* [uj ./• JrL'*\^el 

••iv-.i.f'iyi/ (2, 

-(or - *>^orJy}r 

and the conservation-law form of the convective terms is expressed as 

T 0 ) - l[r( f f ) -’4 ^ ( 3 ) 


The disturbance stream function is defined as Y ~ Y - y where 


2/^'' f 


^(JT V ^ 


jV 


I 






(A) 


The mapping into the unit circle is accomplished by the transfonnation 


where z ■ x + iy 
K ■ re 
Y - 5 + In 
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A proper choice of the conatance y and c Invokaa the aolutlon for Che 
flow over any one of a tivanber of stiapes Including a flat plate, a 
circular or elliptical cylinder, or thick airfoils that may be ayn- 
metrlcal or cambered. The trailing edge of an airfoil shape may be 
rounded off by defining 

0<lil «’> 


In order to minimize the effects of Reynolds number on the uniformity 
of the solution, the radial coordinate is stretched according to 


with 

. - *1 « * _ J ~ )J 

V -■ . K, i I 

The constants r^, k^ and k 2 (all positive) determine the value of p. 

As r varies from r^ to 1, p varies from 0 to 1. The Reyunlds number 

Is R ■ U£/v (i and v are, respectively, the chord and kinematic vlscosclcy) 

and L Is Che dimensionless chord. 

The components of velocity u, and u^ are defined in terms of the 
disturbance scream function as 



,, . I ia-) 


) 


( 8 ) 


For Reynolds nmbers nuch larger than unity, the vurclclcy In the 
flow field exists only naar the body and in the %nke. Away from this 
region, Che flow is essentially irrotaCional. Therefore, Che region 
of calculation is subdivided Into two parts: a small viscous region 

and a large irrotational region bounded by and with ^ 

(see Figure 1). 

The boundary conditions applied along the boundaries of the 
region of calculation are as follows: 

1) On the surface (p ■ 1), the constraint of no slip is applied 
in the form 

and !LY - ^ (10) 

IT r’ 

Condition (10) is used to calculate the surface vorticity from the 
stream function equation (A) with replaced by V. 

2) Along the line of symmetry, the vorticity u and the disturbance 
stream function tfi are specified to be zero. 

3) The flow at the far boundary is constrained with first-order 
differential relations obtained from the liavier-Stokes equations by 
dropping the tangential derivative of the pressure and viscous terms, 
i.e., at the outer boundary, the inertia terms are dominant: 
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( 12 ) 


and U 2 !■ obcained frum the 6 component of the Navler-Stokea aquation. 
That la. 

Note that. In (13), w • 0 wt»en 


(13) 

At t ■ 0, the flow la 
2 


Irrotational (without circulation), l.e., w ■ 0 and ■ -yr . 

Theae boundary conditions are believed to be superior to specifying 
either potential flow or uniform velocity since eddies or vortices can 
pass through the downstream boundary. Also, since the velocity fir 
away Is not defined, the circulation there can change with tine. There- 
fore, equation (11) correctly represents the vortlclty transport through 
the downstream boundary. In equation (13), the absence of the tangential 
pressure derivative will not significantly affect the motion of a vortex 
through the boundary. 

The surface pressure distribution Is obtained by Integrating the 
tangential component of the Navler-Stokes equation. That Is, 


' R jr 



(14) 


where p(0) • 0. The pressure coefficient Is, therefore, equal to 

2p. On the surface, the tangential stress Is given by Oj ^2 “ (L/R)u. 

2 

Both p and 0^2 dimenslonlers with pU . The coefficients of 
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lift, (lr«g, and noniint around Ch« origin of tha z plana (daflnad 
as positive In the counter clockwise direction) ara given by 


IT, 


t 7“ da' 

0^ 


‘ < c -» t. jr j- K 

o c 

( • — * * ; i(_ ^ r 'if 

'u “►* I- I K '■ 


(15) 


(16) 






■M*- 




wnera the subscripts P and F, respectively, represent the contributions 
due to pressure and viscous forces. 

The Newton-Kantorovlch procedure transfoms equation (1) Into a 
sequence of linear partial differential equations (see Apoendlx A for 
developoent) each of which has the fom, 
f(Jr' «J V^» ^ 


r j(^ j-j( '^)j 

^ ~ ‘ io 


where 


(18) 

(19) 


Since the non-linesrlty in equation (1) occurs only In the Jacobian 
term, the Hewton-Kantorovlch transformation affects this term alone. 
The subscript nH-l refers to the variables that must be determined by 
solving the current linear partial differential equation and the sub- 
script m refers to the variables that were determined by solving the 
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previous linear partial differential equation In the sequence. The 
Jacobian tenna on the right hand side of equation (18) are now linear 
becatse only the ''arlables with subscript nrfl are unknown. A sequence 
of these equations are solved such that the solutions and converge 
to the exact solutions and w of the non-linear problem described by 
equations (1) and (4). As shovm In references 1 and 3, the successive 
approximations described by equation (18) converge quadrat Ically. 


t 
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NirMERICAL FORMULATION 


Two conputer programs were used for performing the calculations 
of the current study. The first computer program utilized tlte same 
numerical technique as that used by Mehta^ for solving the non>steady 
incompressible Navier- Stokes equatione Tor arbitrary airfoils at angle 
of attack but was a simpler version only applicable to symmetrical 
slrfoils at zero angle of attack. This computer program is referred 
to as the original or unmodified program and was used as a means of 
checking the validity of the Newton-Kantorovi..h solutions. The second 
progrsm was sl-o*.],';'’ to the first one except that it was modified so as 
to solve t'tM sequence of linear partial differential equations required 
by the Newton-Kantorovich concept and expressed in equation (18). The 
details of the numerical procedure used for the original program may 
be found in Reference 6. The nisnericsl procedure used in utilizing 
tb*'' Newton-Kantorovlch concept is shown here. 

Each Newton-Kantorovich linear systet. was solved by a three-point 
backward time differencing and implicit factored central space dif- 
ferencing scheme. 

The spatially factored, time-differenced, expression for the 
Neifton-Kantorovich version of the vortlcity transport equation (18) 







¥ 


( 20 ) 
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where the operators 6^ and 6^ are given aa 

- i." ♦ £ it £ f -4 1 

k J, ' 3 Jr t y S5 J 

and the croaa derivative terns Q, P, and G as 


- . 1 A , y‘ i ( r 1 1 

VI • - ^ L L^«»l *( ) Jj 

r . J r i ^ K.". •><< ^ 

' * j I L ^ / "^ 5^ 

ja\ n^*> »». / J 

s, xi3£| - ik:±£- 

3 tir L*-*^ c)^r / J>^ J.y J# 

. i u>'’ i^c 

\ (r.», Jc^ ** a / J 

The values of the three-point backward time difference parameters 


T - 3 



and the trauforioation parameter A is 


A 


hV£ 


(21) 

( 22 ) 

(23) 

(24) 

(25) 
are 
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Consistent with the spstlsl factoring concept » equation (20) is 
split into two equations whose finite difference aitilogies each produce 
tri-diagonal systems of equations. That is, 


and 


^ ^ »*i>| 



TM 




w 




f T, w 


r^) 


(26) 


,, if r 

7"^' A A 


(27) 


The spatial derivatives associated %rlth equations (26) and (27) as well 

as equation (19) were approximated by central differencing formulas 

everywhere except at the boundaries of the flow field. The truncation 

2 9 2 

error for the vortlcity equation is 0[(4p) + (A6) + (At) J except 

at the first two time steps %diere the temporal error is 0[At] because a 

two-point backward .'difference fcrmula is required such that T^ “ 2, 

T^ ~ -2, and T^ > 0. The truncation error for the disturbance stream 

2 2 

function equation is 0[(Ap) -f (A6) ]. 

The no-sllp boundary condition represented by equations (9) and (10) 
is reformulated in terms of vortlcity. The finite difference expression 
for this boundary condition is 


''3L 




(a+.3ar)H 


I r(>(^n ) 

■ ■ — I ■■ I m^i 

■‘L 





A ^ 

Ar iJr 


-3 





(28) 
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The truncation error for this formulation la 0[Ar ]• In addition, the 
condition t ■ 0 la alao aatlafled on thla lioundary. Along the llnea 
of symmetry the boundary conditions u ■ 0 and i(i ■ 0 are apecified. At 
tho outer boundary, the finite difference forma of equations (26) and 
(27) are used to describe the vortlcity condition except that the 
diffusion terms are deleted and forward space differencing is used 
in the p direction. This fornulation results in a first order trun- 
cation error in the space variable p. The disturbance stream function 
at this boundary is obtained from the finite-difference form of equa- 
tions (12) and (13): 

rr,»/ rr,,i J 

1 (30) 

The present values of the vortlcity at a wall grid point is determined 


tAj” +£ (CiUj) 

m.i m* \ II \ / 


where is a relaxation parameter and k is an iteration counter. 
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Once the solution to the linear syeten of aquations has been 
obtained as described above at a particular time step, the Newton- 
Ksntorovich procedure is invoked to produce a new linear systen and 
the process is repeated as many tines as required for convergence to 
the solution of the non-linear system for that time step. Then the 
entire procedure is advanced to the next time step. The initial guess, 
w^, at the solution fot each time step, was determined on the basis of 
an extrapolation of the converged solutions at the two previous time 
steps. That is. 


CJ 


• ►h' 
0 




r> 

U) - 




(32) 
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This approximation has a truncation error of 0(At ). At the first 
time step was set to zero and the potential flow stream function 
was prescribed for At the second time step and i|/^ are set equal 

to the converged solutions at the first time step. For subsequent 
time steps, t|i^ was set equal to the converged solution at the previous 
time step. 

Since the truncation error of the Newton-Kantorovich approximation 
is 0[Au • Aij)], the truncation error of the initial guess at the first time 

step is not determinable in terms of At, but at the second time step it 

3 4 

is 0[At ] and at subsequent time steps it is 0[At ]. 

Computations of pressure coefficients on the surface and a deter- 
mination of the loads are made for each time t with a finite-difference 
integration fornula derived by cooibinlng two four-point expressions 
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(In order Co have a lower effective truncation error). The integral 
of any function f (p) between the nodal pointe 1 and 1 -f 1 on Che 
surface are expressed as 




I + 1^ 



where j ^ J. '♦’*♦ 3 ) 

This formula is used first to determine the pressure from equation (14) 
and Chen is used to compute the loads. The pressure calculations 
require the normal vorCicity gradient, which is represented by 



Tills equation has a truncation error 
truncation error is not used because 
vorticity with a truncation error of 
identical to the above equation for 


of 0 (Ap). A formula with a smaller 

Che derivation of the surface 
2 

0(Ap) requires the use of a formula 


I 


I 


I 
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DISCUSSION (f NUMKKICAl. RESULTS 

The two computer programs discussed In the previous section were 
both used to determine the flow field about a circular cylinder at a 
Reynold? number of 15, and a 12% thick symmetrical airfoil at zero 

3 4 5 

angle of attack at Reynolds numbers of 10 , 10 , and 10 . In all 

cases, the bodies under consideration were impulsively started. The • 

cases for the circular cylinder at a Reynolds number of 15 and the 

airfoil at a Reynolds number of 10,000 were computed until a steady 

state was achieved, while the other cases were run for early time only. 

*^6 results calculated by both programs were compared with each other 
as a basis for determining the number of Iterations of the l^ewton- 
Kantorovich procedure required for convergence. Figures 2 through 5 
Illustrate some of the numerical results and table I suomarlzes the 
convergence characteristics of the technique for these calculations. 

In general, the results obtained by both techniques were in excellent 
agreement with each other. As seen in the figures, the maximum dif- 
ferences between the two sets of results was about 1%. This difference 
did not decrease with an increase in the number of Newton-Kantorovi h 
iterations but did decrease with a finer spatial grid. This indicates 
that the difference between the solutions was due to the slightly 
different finite difference approximations to the Navier-Stokes equa- 
tions produced by the two techniques. 
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In addition. Table I Hhowa that finer apatial gride were required 
at the higher Reynolda nuvbera in order to achieve the Kame level 
of agreement. This was due to the fact tl d the aolutiona berome more 
non-uniform in the circumferential direction with inctcaaing Reynolds 
number. Since the ciremferential coordinate was not stretched %rlth 
Reynolds number ns was the radial coordinate, the number of grid elements 
in the circumferential direction Ivad to be increased in order to maintain 
the same level of accuracy in the regions of large gradients in the 
solution. 

The convergence characteristics of the Newton-Kantorovich procedure 
were excellent for all test cases performed. In fact, after the first 
several time steps of each test case, the initial guesses for succeeding 
time steps were sufficiently accurate that only a single iteration was 
required to achieve convergence; and thus the technique became completely 
non-iterative with respect to the non-linearity in the equations. This is 

4 

illustrated in Figure 6 for the Reynolds number 10 wliere the vorticity 

at a particular point in the flow field is plotted versus time using 

one Newton-Kantorovich Iteration and is compared to the same case using 

two iterations. For the first several time steps, two iterations were 

3 4 

required for convergence of the Reynolds number 15, 10 , and 10 cases, 
and a third iteration was required for convergence of the Reynolds 
nimiber 10^ case. These succeeding iterations at the earliest time 
steps were required because the initial guesses at the solution for 
early f:lme were very inaccurate (see Nunerical Formulation). 
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Smaller time Btepu were required at the higher Reynolds numbers 
to achieve a sufficiently accurate Initial guess at the solution to 
preserve the number of iterations required for convergence of the 
procedure. However, smaller time steps were required anyway to achieve 
convergence of the surface vorticlty boundary condition at these 
Reynolds numbers. Hence, these time step restrictions are consistent 
with the boluticr procedure and not a limitation due to the Newton- 
Kantorovich process. 
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CONCLUSIONS 


Basod on the results obtained In the current study, the Newton- 
Kantorovich technique can be succcsafully applied to the nuiaerica] 
solution of the laminar, non-ateady, incompressible, two-dimensional 
Navier-Stokes equations. The convergence characteristics of the tech- 
nique were excellent for all geometries and Reynolds numbers tested. 

In fact, except for the first several time steps, the procedure requires 
only one iteration to achieve suitable convergence; that is, the pro- 
cedure becomes non-iterative with regard to the handling of the non- 
linear terms in the Navier-Stokes equations. This conclusion indicates 
a potential for significant reduction in computation time over other 
current iterative techniques. 
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AHPtNDlX 

The nonlinear, mmatatlonf r> Navler-Stokeh oquatlonN, which can 
be described niattieinaclcally by 


<^{^3] = I , 1 , ^ 



with 


( 1 ) 


can, equivalently, be described in the form 

P(^)=o 

/J. 

Operator P of (2) will be considered a mapping from the Banach Space C^(S 
into the Banach space C^°(S)» where C^*^(S) ^(x,y,t)js times continuously 
differentiable in (x,y)€ S,T times continuously differentiable in 
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t,0^ c <^, Inciuulve of special derivatives up through order 


M} 


Not a 


ti^t ( Jdenotcs greatest Integer part. (S) will be considered under 


the norm 


ll'fll - 


C 

T 


tz » 


l¥j 




I ^ 





a- 



In order to apply Newton's method ns generalized to function spaces 
by Kantorovich^, It Is necessary to determine whether or not P', the 
Frechet derivative of operator P of equation (2), P", and (P') ^ exist, 
and if they exist, upper bounds for/^P//, /(p'y/, and|^(P') ^/|/ imist be 
determined. 

It can be readily shown (see (3), (4)) that P'(>jj) exists for 
arbitrary Cj^(S), and that 







y/S'Vo 


Similarly, it can be shown ((3) (4)) that 

for arbitrary lii^, 0 in Cj^ (S) . Also, 

“ H f’ix)ll H H^// P'Y't) J! 

4 

for some initial guess ^ (S) t then the algorithm 




I 


24 




r*o 


convcrgui to the solution of (2) (>.se (5)). 


-I 


Thsrefore. It retnains to detsnnlne If and wiisn P'(v ) exists and 
* o 

to obtain an upper estimate for jj ** Jf 

To this end, let m ■ 0 In Eq. (3). Therefore, 

F(X) -h c: 

Let ^ ; therefore, - If 0*' exists, 

- F(^) p{x) . 

Hence, we seek a solution^ of the problem 

o = PCX) -*pYx)T, :,o_W 

t^o 

Let PV'^} = Y/4 ^ 

(oiftee A ij ^!<*€Ar C.a (S ^ , 

An)^ = 4_ 

Equation (4) can, therafora, ba axpreaaad by ( 4 ^ 


43'*F - 


*'AA y- an - 


e^X)- 4 lx )x, aVAOc.^ 


si (5) 

o 


J.(i) (l) 


then 


25 


/ 


Lat ^ b« defined ae followat 


^ - (H V;»;v>’^ 


where 


It can be directly shown that It ^ of (6) exists, then ^ satisfies 
(4), where G Is the Green's Function of La Place's equation for S 

y 

and H the kernel function for the heat equation 

^ <p =■ o , 5 


(2<^)=o ^4 ~ ^ 

r»o 


In order to determine sufficient conditions under which ^ c'f 
equation (6) will exist, let B denote the linear operation: 

Equation (6) therefore takes the form 

(x- qa)^ = $ ph() 

3 3 

Let Cq (S) S the nompletion of Cq (S) under the norm 


2h 


M hi t hi 


o* t 


^ajf 


<Xit<-^ 


V 

a* r- 


l%l 


V-i 

o% f <*d 


;,0 ^ 


V 

iJfc f <.*A 




-f If^A t I ^'V'l I 


IK* l^'^l , 


• t-t <.^ 


o» f <aO 

^3 ^1 

If (HA) is considered a pupping from Cq (S) into C^ (S), and 

jl BAj^4l, then (I - JiA) ^ exists, and^ • (I - BA)”^BI’(^). Therefore, In 

order to show that 'i* exists. It Is sufficient to show that llelilMiUI 

^ i 0 /u ^ r® 

To this end consider A as a mapping from C ' into C . For i ^ ^ , 

" 0 o / 


h'i'Jl ~ ‘K«ur/'V’/ \ f«r M ^ C>' 


U 




/« 


; (I'ThMilXl Wv.JtlX 1 

Ole*4>»s 


Since //All “ ^ 
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by definition. 


A/'VO </ r ^ ^ ^4>-+ 4% 4,j 
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1 


t 


I 


21 


and 


II ^ ^ |-''^ ^4jf '~ ^ I; 


It follows that // ^If^Al , wtiere 

X 


/M = i/>v^ a ( I'^^lj 1‘^oJj / ^'% \^ 1^"%^ ) ■ 


X V‘S 

o i -t ^-#6 

o 3 

Consider B as a mapping from C Into Cq . 


by definition. Tlierefore, since 


I Its II- Vf (i}f^ IL \ 

. I _ ii vl'* ' 


Mr 


IISi-lL- 1 1 Kf»k/f| h#)^ / 

I ( 

^ ^ ^ ^ ^ Ul^ I 


/ 


OLr«.»4 

-♦• 


-• /6k^n<v.“^*/4s<^' 

rj ^ /'^r% v> -^V v: ' 

Le-t i 

r. JKA» Ti j" /^ V ■■, 1 r){. <• +/«/r/^/v /■* /^ /»/<^/ 

J 
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Tlierefore, for ^ cxistu, and ^ ■ (I - BA) BH(^). Hence, 

It follows directly that -P*( q) ^ ^ that 

y //^ll /t4 

Froci the preceding, it directly follows that II^’q, the Initial 
educated gueas. Is judiciously chosen, a sequence of functions 

■ ], 2, . . .) can be directly constructed that conver^ the 
exact solution of the nonlinear, nonstationary Navler-Stokes problem 
as expressed by Equation (1). In addition, an explicit error estimate 
can be directly determined for each approximate solution . 


Explicitly then, for 


II ] // 






with , the sequence of determined by Equation (3) exists 

and converges to the exact solution 'V** of Equation (1). 

Moreover, under these conditions, the method yields the explicit 
error estimate 
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The actual application of the Newton-Kantorovich method to equation (1) 
produces a sequence of equations of the following form: 


'i / 

c » ' 

Define the vortlclty, u, as 

UJ^ - 


Hence, the final form of the equation becomes 


*j( iif) 

Transforming thla equation Into the interior of the unit circle and 
transforming the radial coordinate to a new variable p that allows for 
stretching of the solution in that direction produces 


H 7c ‘ /if 
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Table 1 Summary ot test caseH performed durint; Che current study. 






4 



Fig. 1 Domain of calc ition, boundary conditions, and grid notation 
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Surface vortlcity and pressure coefficient distributions on a 
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Surface vorticity and pressure coefficient distributions on a 
12Z thick symmetrical airfoil; a ■ 0“ : Re ■ 1.0 x 10^. 
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